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ABSTRACT 


The non-linear elasto-plastic response of a clamped 
flat plate to a typical air shock wave were investigated. 

The nonlinear effects to the plate responses due to the 
material and geometric nonlinearity were studied. 

In this study, (1) the necessity of the modification of 
old armored vehicles was reviewed and (2) the NASTRAN code 
was employed in this investigation, (3) the theoretical 
background of the nonlinear transient analysis was described, 
(4) a step by step procedure of analyzing the dynamic load 
problem by shock wave using NASTRAN code was developed, and 
(5) sensitivity analyses were performed and also difficulties 


associated with the nonlinear analysis were described. 
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ft. glitWeRODUCTHAON 


Conventional tanks (old tanks) which have steel armor, 
can be completely defeated by one round of modern projectiles 
such as anti-tank missiles. That has been proved in the 1973 
Arab-Israeli war. In that war, the 40% of total Israeli 
casualties were caused by armored vehicle destructions. Thus, 
tank armor improvement became the first priority in the Israeli 
tank corps. 

If war occurs in any country, the same problem will be 
faced. That is the reason why the modification of the conven- 
tional tank armor is needed against modern anti-tank weapons. 
To improve the armor protection, the material properties, 
projectile characteristics, and explosion phenomena must be 
understood well. This is a complicated and difficult problem 
and is originally a nonlinear problem. 

First the shock wave from an explosion imposes a dynamic 
load on any object. That dynamic load is characterized by 
rapidly reached peak overpressure which then decreases as the 
shock wave decays. The net effects of the load depend both 
On the nature of the shock wave and the geometry and construc- 
tion of the object. For example, if the explosion occurs on 


thin skirting plate, perforation takes place by impact load 


and their impulse loading for main armor would be reduced 


significantly. The typical air shock wave was characterized 





as follows; the shock wave 1s supersonic and generates the 
abrupt increasing overpressure. 

The non-linear elasto-plastic responses of a clamped flat 
plate to a typical air shock wave were investigated. The 
nonlinear effects to the plate responses due to the material 
and geometric nonlinearity were studied. 

In this study, (1) the necessity of the modification of 
old armored vehicles was reviewed, (2) the NASTRAN code was 
employed, (3) the theoretical background of the nonlinear 
transient analysis was described, (4) a step by step proce- 
dure of analyzing the dynamic load problem by shock wave 
using NASTRAN code was developed, and (5) sensitivity analy- 
ses were performed and also difficulties associated with the 


nonlinear analysis were described. 





Tl. LITERATURE REVIEW 


A. MODERN PROJECTILES 

An important version of the AP (Armor-Piercing) round is 
that known as APDS (Armour Piercing Discarding Sabot). [In 
this round the armour-piercing element is a hard core of 
calibre significantly smaller than that of the complete 
projectile and shaped for maximum anti-armour effectiveness. 
This is surrounded by a much lighter casing, or sabot, which 
breaks up as the projectile leaves the gun barrel. The 
point of this arrangement is that the core can have a shape 
which differs considerably from that of the complete projec- 
tile and, not having been scored by the rifling in the barrel, 
is aerodynamically clean. Moreover, since it is very much 
heavier than the sabot, the core acquire a greater energy 
than could be imparted to it if it were fired, without sabot, 
by a gun of appropriate calibre and the same chamber pressure 
as that of the larger gun. A controversial refinement 
involves the use of depleted uranium in the core. This 
metal has an exothermic reaction with steel which aids 
penetration; it is controversial because of its connection 
with the other isotopes of uranium and their radioactive 
properties. 

High explosive (HE) rounds with contact fuzes can also 


be used against tanks and may, if they are large enough, be 


Ie, 





effective against exposed portions of the tank mechanism 
(such as the tracks) but will not in general be effective 
against heavy armour. Two modern refinements of the HE 
round may be effectively used in some circumstances: they 
are the HEAP (High-Explosive Armour Piercing) and HESH 
(High-Explosive Squash-Head) rounds. 

The HEAP round has a projectile of relatively high 
velocity, a hardened case and a high-explosive charge which 
1s contact fuzed with a very short delay. The effect of the 
arrangement is that the projectile either embeds itself in 
the armour or passes through it and then explodes with an 
effect which is more damaging than that of the simple HE 
round. 

HESH rounds can be used with lower-velocity weapons, and 
are so designed that the nose of the shell crumples on impact 
to make a large contact area before the explosive is detonated. 
Such rounds will not, in general, pierce armour, but the 
explosion sets up shock waves in the armour which in them- 
selves can incapacitate the crew and may also cause pieces of 
metal to break away from the surface remote from the explosion. 
Such fragments fly off at high velocity and can do great 
damage inside a tank turret. 

Most widely used in infantry weapons, however, is the 
projectile (which may be a shell, a rocket, a grenade or a 
guided missile) with a shaped-charge warhead frequently 


Gescribed as a HEAT (High Explosive Anti Tank) round. This 


aE 





type of projectile is designed to be effective with only a 
moderate impact velocity and depends for its effect on the 
design of the warhead and is formed round a hollow cone with 
its apex to the rear. The base of the cone is of approxi- 
mately the same diameter as the warhead at that point and in 
front of it is an air space enclosed by the rounded or 

pointed head of the projectile on which the fuze igniter is 
mounted. The operation of the fuze is critical to the success 
of the charge. It has to have a set distance from the surface 
of the armour when it is detonated and for this reason the 
tae ter 1S Carried well in front, often giving rise to the 


distinctive pointed nose of a hollow charge round. 
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Figure 2.1. The Effects of Current Anti-Tank Projectiles 
on the Steel. 
The igniter is connected to the detonator which is at the 
back of the charge, this starts a detonation wave running 
meameuga the conical explosive and the result is that the 


force of the explosion is concentrated into a thin jet 


eZ 








fering wat dad aspeed Of Over 20,000 m/s, and at a great tempe- 
rature. The jet takes a finite distance to form, hence the 
importance of having the igniter in front and the distance 
required to form a jet is known as the "stand off". For most 
modern warheads the stand off distance needs to be between 2 
and 3 diameters of the warhead and to provide the full dis- 
tance sometimes interferes with the ballistics of the 
projectile. 

The effects of a hollow charge on the interior of a 
vehicle can be variable. A very small jet may well not 
Cause much damage apart from making a hole in the armour and 
destroying whatever 1S in its direct path, but it is possible 
to “turn” warheads so that the penetrative effect is lessened 
but there is more energy left after passing through the 
armour and this energy can be spread out into more of a fan- 
Shaped exit blast with appropriately enhanced damage. All 
hollow charge jets carry a certain amount of molten metal 
into the interior and some research has gone into finding 
ways of carrying more through the hole. 

Two examples of the Soviet anti-tank guided missiles are 
described below. [Ref. 1] 

ie oeogesmamei-Tank Guided Missile 

Sagger is the NATO code-name for the Soviet anti- 

tank missile which is also known by the US alphanumeric 
designation AT-3, which is now extensively deployed in 
Europe and elsewhere. Currently it is known to be used in 


at least three ways: 


of 





a) 


b) 


Cc) 


As a man-portable missile 
Mounted on a bracket 


Mounted on a vehicle 


The characteristics of sagger missile is as follow. 


a) 


i) 


Type: Wire-guided MCLOS (Manual Command to Line 
OLesaene) 

Welgmt at launch: Il1.3 kg 

Diameter: 120 mm 

Warhead: HEAT 

Range:  200=3000 5m 


Penetration: 400 mm 


2. AT=-4 Anti-Tank Guided Missile 


There have been reports for some years of a Soviet 


SACLOS (Semi-Automatic Command to Line of Sight) system, but 


Dees NOt until 1980 that photographs and a little informa- 


tion appeared in some Warsaw Pact military journals. 


Siaracteristics: 
a) Type: Wire guided SACLOS system 
By Werght at launch: 40.4.kg 
c) Diameter: 120 mm 
d) Warhead: HEAT 
e) Range: 2000 - 2500 m 
f) Penetration: 500-600 mm 
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Figure 2.2. Sectioned View of Sagger Missile. 


Which of these various types of projectile will be most 
effective against armour which depends on the design and 
thickness of the armour. Against solid armour the HEAT 
round is particularly effective; but in certain circumstances 
a sandwich construction is likely to be more vulnerable to an 
APDS round. There is therefore a good deal to be said for 
having a mix of projectiles available for satisfactory 


defense against tanks. 


B. CONFIGURATION OF ARMORED VEHICLES 
1. Armor 
We can define the military armor as a material 


intended to protect men and installations of equipment. 


iS 





Aluminum, titanium and magnesium alloy, plastic laminates, 
nylon fabrics and masticstone composites have all been used, 
but the most common form of armor is heat-treated alloy steel. 
This comes mainly in two forms: homogeneous and face-hardened. 

The first has constant ballistic and mechanical quali- 
ties throughout, which in modern tank armor means high tensile 
strength (of the order of 150,000 lb sq in) and medium hard- 
ness (i.e. Brinnell N 300); face-hardened armor consists of 
a very hard first layer which resists penetration but lacks 
toughness and a more elastic core to support the face. Naval 
deck armor and modern tank armor are homogeneous, while face- 
hardened armor is still used for some light armored vehicles 
and side protection for naval vessels. 

Steel armour usually means rolled or wrought plates, 
but castings are used for tank turrets and even for entire 
tank hulls, as in the M.48 and M.60 tanks (Steel armor is 
usually a nickel-chrome-molybdenum alloy, but other alloys 
including manganese~-silicon-chromium-nickel steel have been 
used, especially in Soviet tank armor.). 

Typically, a well-designed hollow-charge warhead can 
be expected to penetrate armor up to five times the diameter 
of the warhead. The best defense against such attack is not 
to make thicker armor, but to dissipate the jet before it 
reaches the surface. To do this a light screen is carried 
some distance in front of the main armor, of sufficient 


thickness to set off the fuze of any hollow-charge projectile. 
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OuLee naturally such screens do mot last long, but they are 
easily replaced and provide almost complete protection for 
the armor behind. A common use of such a technique is the 
thin skirting plates carried over the tracks of all modern 
tanks. Before leaving this subject it is appropriate to 


mention a recent British armor development, the 'Chobham' 
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Figure 2.3. The Effect of Anti-Tank Projectiles on "Cobham" 
Armor. 
armor. This is claimed to offer greatly enhanced resistance 
to all forms of projectile as is indicated by the accompanying 
illustrations. Details of the construction of this armor have 
not been released, but it seems probable that it is a sand- 
wich made up of two sheets of armor plate separated by a 
highly dispersive, possibly granular, filling. The disad- 
vantage of this armor lies in its bulk. To gain the protec- 
tion, it has to be thick and fairly heavy; tanks are already 
large enough and the addition of armor plating, which is 200 
mm or more in thickness makes a startling difference to the 


outline and bulk of the entire vehicle. 
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2. Curment Medification of Old Tanks 

Determining which one should we choose, a new one 
or modification of an old one usually depends on cost and 
Powatical situations. 

The U.S. completed conversion of M48 tanks to M48A5 
standards in 1980. This conversion basically included adding 
diesel engines to M48Als and adding 105-mm guns to the tur- 
rets of M48Als and M48A3s. Numerous minor changes were also 
made. But little improvement in armor protection has been 
made. If we want to take more modifications for M48A3s and 
M48A5s tanks, a spaced applique armor should be added and 
include a 1,200 hp engine, improved transmission, better 
fire control, and a hydropneumatic suspension system. 

The following example shows modern modifications. 

The High Performance M60 MBT has been developed as a 
private venture by the General Products Division of Teledyne 
Continental Motors. It is essentially an M60 series MBT 
Fitted with a new powerpack, new suspension and, as an option, 
additional armor. These modifications can be carried out on 
existing M60 series. Teledyne believes that to fit an M60Al1 
with the new powerpack would cost under $200,000, while the 


added armor would cost just $40,000. 
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Figure 2.4. Applique Armor in M60 MBT. 


3. Israel Add-On Armor 

Israeli experience in the 1967 campaign proved that 
mobility was not substitute for armor protection and they 
therefore decided at an early stage that the main emphasis 
would be placed on armor, with firepower and mobility second 
and third priorities. After that they spent a lot of money 
Beeeersure OUL this protection problem. Total cost of 
research, development, trials and the building of the proto- 
type vehicles was about $65 million. 

In May 1977 Israel finally announced that it had 


developed a new MBT called the Merkava which was made of 


JES 





cast and welded armor with a well-shaped glacis plate. 
Behind the first layer of cast armor is a space filled 
with diesel fuel and then another layer of armor. This 
Spaced armor gives the tank protection from HEAT projec- 
tiles and ATGWs (anti-tank guided weapons). 

They also modified their old tanks such as M60s, 
Centurions, etc. 

The add-on and readily detachable protective arrays 
seen on Israeli M60s and Centurions contain ceramic tiles, 
would also mark a first in terms of operational deployment. 
The possibility cannot be excluded that armor of both kinds 
1s involved. 

While no answers are to be found in external appear- 
ance, the configurations, as seen in the accompanying Defence 
Attache photographs, show obvious differences in the thick- 
ness and shape of the packs provided for M60 and Centurion 
respectively. The Centurion's additional armor is mainly in 
the form of quite shallow panels fitted to the brow of the 
turret and thicker pieces to the glacis and the mantlet 
either side of the main gun; a hole is included in the left- 
hand piece to permit the firing of the coaxial machine gun. 
Wedge-shaped packs are mounted on top of the storage boxes 


fitted forward above the track guards. 
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Mure 2.5.. Add on Armor Blocks in M60. 
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Figure 2.6. Addon Armor Blocks in Centurion. 


The M60, partly because it lacks whatever protective 
values come with the Centurion's storage bins, has been 


provided with a considerable amount of add-on armor around 


22. 
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Eaee turret... Here the packs ane thick.and vertically deep, 
with smaller and thinner sections patched in -- with no 
great concern for aesthetics -- to the upper surfaces which, 
because of the turret's curvature and varying angles to the 
perpendicular, not to mention obstructions such as lifting 
eyes, cannot be covered by the main boxes. As with the 
Centurion, panels are fitted to the front of the turret roof 
and to the glacis; in the latter case, because the M60's 
basic glacis is curved, it has been necessary to fita 
straightening frame to which the packs can be attached. The 
M60's installation is completed by the fitting of wedge- 
shaped or angled sections either side of the driver's hatch 
to help protect the vulnerable area where turret meets hull, 
with further angled sections each side of, and thinner plates 
above, the mantlet. 

Although it can be stated with a fair degree of 
Semeaioty that the add-on arrays are imtended to help defeat 
shaped-charge high-explosive anti-tank (HEAT) attack, any 
discussion of what the packs and panels consist of must, in 
the absence of any word of guidance from the Israelis - who 
on the contrary are happily watching, if not promoting, the 
Circulation of a number of often fanciful and frequently 
fimetally exclusive ‘solutions’ to the puzzle-fall short of 
a definitive answer. 

While one source suggested that the boxes were empty 


(and the large ones on the M60 do sound quite hollow when 


ZS 





haumped) the: probabilaty us that there is» mere to it than 
that. In view of solutions that have been adopted else- 
Where, though usually in the context of new-build tanks, it 
seems quite likely that the Israelis have developed their own 
Gemeamic tiles» and that the boxes are partially filled with 
them. The alternative of a form of active armor is suggested 
more by the vehemence with which several sources in Israel 
stated this to be the case than by any objective evidence. 

A number of different types and configurations of 
active armor are possible. An example is a sheet of explosive- 
laid of top of a steel plate, with a covering for environmental 
protection on top - which detonates when hit by the incipient 
jet of a shaped-charge warhead and in the process disrupts the 
jet's further formation. If the Israelis are indeed deploying 
a form of active armor, such a fact would carry the implica- 
tion that they have sufficiently overcome the formidable 
difficulties that seem to have prevented larger and better- 
endowed countries from doing likewise; the United States and 
Soviet Union are among the nations that are known to have been 


engaged in related research and experimentation for many years. 


C. ELEMENTARY DETONATION THEORY 
1. The following assumption will be taken for the 


simplest analysis of a detonation: [Ref. 2 ] 
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S|) Gren i Low “netlonwofimthe shock front) is one 


dimensional 

b) The plane detonation front is a jump discontinuity 
at which chemical reaction 1s assumed to be 
instantaneous. 

c) The discontinuity 1s time-independent, so that the 
product composition (at least directly behind the 
shock front) is likewise time-independent. 

u represents the particle velocity of blast wind, which is 
@iemeor the products behind the shock front. Subscript x 
represents the initial values of the properties, y those 
after shock. D is the detonation speed and p is the shock 
Wave overpressure. 


From mass conservation across the front we can write 


D= D- i 
O = ( u,) (1) 





Figure 2.7. Explosive Shock Front. 
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and from conservation of momentum, 


Py Py #9, 4,0 (2) 


Eliminating up from the two equations we obtain the equation 


for the Rayleigh line: 


Ze 
OD” - (py - P/ (vy > Vy) 


ll 
Oo 


(3) 


(Observe that in part of the previous equation we replace 1/0 
by the specific volume, v). Since for a detonation, ae 2 
it follows that Sy a (for real values of the detonation 
velocity); and, indeed, it is observed Bea rule that there 
1s about a one-third increase in density on detonation. If 
MreeeecOudition 1s substituted in eq (1), it follows that 
OF < u, <d, so that the products travel toward the detonation 
meome, though trailing behind. 

If the Rayleigh equation is solved for Ey we obtain 

Qe 


— Z it 
Ee = (Dy + PD ) Ox “D a (4) 


wherein op 1s linear with Yoge the slope of the line being 


-— aD. . 


- The conservation conditions thus require that the 


initial and final states of a detonation must lie at the 
termini of the Rayleigh line. 
Conservation of energy leads to 


= 7 ae _ 2 
ey’ eo ee > Us = E(P,1/x,A=0) +P vv t+1/2D (5) 


in which i = the extent of chemical reaction (A =0 for no 
reaction and 1 for complete reaction). Eliminating wu _ and 
Perc lomeGba (5), Saing Eqs (1) and (2) leads to the Hugoniot 


Curve: 
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ey >! = (Pl svi, A=0) _ hyo (PtP) ing a (6) 
AE, = 1/2 See ee (7) 
AE, 1s the internal energy of the detonation products at the 


very high pressure and temperature of the detonation state 
less the internal energy of the unreacted explosive at ambient 
temperature and pressure. It may be written as the sum of two 


terms: 
AE, = AED + AE, (8) 


where AE, is the internal energy change for the conversion 





of reactants to products isothermally, and aES 1s that 
required to warm and compress the products from ambient p 
and T to the detonation values. AE Wn is often symbolized 
as -Q in explosive literature (Q being the heat evolved in an 
explosion when corrected to isothermal conditions). If the 
Pagemcts were ideal gases, then AE. would be given eran, 
and this approximation is often made, though the products at 
pressures of hundreds of bars, are far from being ideal gases. 
2. The Hugoniot curve is that of a hyperbola in the p-v 
plane. For the non-reacting case (A=0) the hyperbola passes 
through the initial state, PyrVyi but for the case we are 
considering here (A=1) it does not. The conservation laws 
require that the final state, however, lies on the Hugoniot 


curve as well as on the Rayleigh line. We consider below 


several cases: 
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For a too small D(=d.) the Rayleigh line does not 


intersect the Hugoniot and thus no solution exists. For 





Figure 2.8. The Hugoniot Curve. 


large D (=D)), there are two intersections, called strong 
(labelled S) and weak (W). For one particular value of 
D es) the Rayleigh line is tangent to the Hugoniot. 
It can be shown that the relative particle velocity, 
a! 1s subsonic at §S caus < ay) Supersonic at W and 
sonic at the point of tangency, called the Chapman-Jouguet 
DOrnit : So oe = Bo (sonic velocity in the shocked medium). 
There exist a number of computer codes based on such 
non-ideal gas treatments (RUBY and TIGER are two of these) and 


moderate success has been achieved in computing detonation 


properties from these codes. The equations are: 
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Py = Kp, p (9) 
ye owen Bo.) (10) 
} 
> = NM%Q? (11) 
Mer 80.762, A =| 22.3, B = 0.0013 (SI units). 0, the 


heat of explosion in J/kg; N, the number of moles of gas per 
kg; and M, the average molar mass of the gas (kg/mole are all 
found from an assumed stoichiometry of decomposition, combined 
with thermochemical principles. Although N, M, and Q are all 
strongly affected by the assumption made as to stoichiometry, 
(femeurns Out that the important parameter, $¢, 1s rather 


insensitive to the assumption. 
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TIT. ANALYSIS OF DYNAMIC RESPONSE OF PLATE BY SHOCK WAVES 


A. THEORETICAL BACKGROUND 
E- eienlinear Static Analysis 
The fundamental method used to obtain a static solu- 
tion is to minimize the error vector {6}, given by the 


following equation: [Ref. 7] 
{6} - {P}+ {Q}-{F} G2) 


where {6} is the error vector of unbalanced forces acting 

at all grid points, {P} is the vector of applied external 
loads, which may change with displacements, {Q} is the unknown 
vector of constraint forces due to single and multipoint con- 
straints, and {F} is the vector of grid point forces generated 
by element motion and stress. The terms are functions of 
displacement, temperature, and stress history. 

It must be noted that degrees of freedom not involved 
in constraints produce null terms in {Q} and that dependent 
constraint points produce no errors. Therefore, when Eq (12) 
1s reduced to the solution coordinates the constraint forces, 
{Q}, disappear. Terms in the vector {F} are dependent on the 
deformations of the finite elements. In linear static analysis 
the vector {F} becomes 


eee “7 [KP ]{u} = {p*(T) } abs), 
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where [kP ] = The linear stiffness Watrix 


{P=} = The "thermal load" vector, i.e., a 
vector of grid point which produces the 
Same element strains as the temperature 
CasSteremeion, 

{uj} = The vector of grid pownt dusplacements. 


To obtain a nonlinear solution, a "tangent" matrix 


{K] is used in which the terms are derivatives in the form 


OF. 
K.. = 


i5 iso, fu =u ale 
J 


1G 


The matrix [K] contains both geometric and material non- 
linear effects. The approximation occurs because the forces 
are highly variable with displacements and nonsymmetric terms 
are approximated. If {ul} 1S a known displacement which 
produces an error {é6(ut)}, then Newton's method may be used 
to predict a new vector (oe! with a smaller error. For 
a displacement {u} near a known solution, fur}, the non- 


linear force is approximately 


{P(u}} {F(ul) } 4 imitasul} ‘Gus 


The substitution of Eqn (15) into Eqn (12) with {6} = 0 
yields the result 

[K] {u-u>} ~ {P} + {Q} - (F(u‘)} fo) 
which can be written in terms of the error {6} to provide 
the format of the Newton-Raphson iteration method: 


Pace = (6 (u")} (17) 
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When [K] may be inverted, Eqn (17) may be used to solve 
fox nes = {(u)}, a new estimate of the solution. 
Equation (17) provides an incremental form for solution 
iteration. Another form of the iteration equation is 
Obtained from Eqns (12) and (17): 


[K] {u***-u*} = {p, we Pee) = fe} = (KP tu} (18) 
where 
feo er teat) } = [R] {ut} + fp} = Ap(a2) (19) 


[K+] is the nonlinear "reference" matrix, available from the 
elements and Ap (u*) 1s the change in loads due to grid 
point motions. The advantage to this form is that the second 
Order "corrective" load vector {£+} equals zero when the 
forces {F} are linear and when [K] contains linear elastic 
terms (see Eqn 13). Because the {f} forces are calculated 
at the element level, all known linear elements and super- 
elements may be bypassed in their calculation. This is very 
efficient when a major portion of the structure remains linear. 
However, the advantages for using Eqn (17) are that the load 
error and incremental displacements (hie ut} are directly 
available for testing convergence. Therefore, in MSC/NASTRAN, 
both Eqns (16) and (17) are used internally. 

In order to avoid calculating the vector {F(ut)}, 
Eqn (19) is solved for {F (ut) } and substituted into Eqn (12) 


to obtain 


Pee ein + poh] {f°} — [K*] {ut} (20) 


oe 








Subtracting ete for two successive values of i in 
em (20) to elaminate the constants a simple equation is ob- 
tained for the load error in the form 

emt ie = eS} - OE) fue - uw 4} (21) 

This form provides a simple iteration procedure, namely 

a) For the first iteration after a change in the load vector, 
Wee Eon (17) and (20). Note that the initial values for 
starting from rest are {f°} = {u°} = 0: calculate 

{§°} = Berens {pt} and solve ([K] ee} Seles (22) (23) 
b) then compute gt directly from tne element and load 
tables: 

- 


eee fk) furl} + fo '} = fAp(ut)} (24) 


c) then compute oe Peem Bem, 21). 


pees (673+ 6°} = kt pata? - u°) (25) 
where {f°} = 0 by assumption. 
d) Use Eqn (17) with i=l again to compute {u-} 


[K]{u - ea = ste (26) 
e) Repeat steps b, c and d to convergence, incrementing the 
Superscripts by one at each iteration loop. 

2. Nonlinear Transient Analysis Using MSC/NASTRAN 

SOL 99 which is a rigid format of MSC/NASTRAN code 
for nonlinear transient problem, is the dynamic comparison to 
nonlinear static systems, providing for combined material and 
geometric nonlinear analysis in the time domain. SOL 66 which 
is a rigid format for nonlinear static problem and SOL 99 


share the same element computer code and solution techniques 
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ama provide Similar data storage and restart facilities. 

The transient solution is performed ina stepwise manner with 
time steps replacing load steps and with effects of mass, 
damping and unsymmetric matrices added to the stiffness 
matrices. 

Because the calculations are identical, many element 
and material discussions in the previous sections also apply 
to the transient analysis system. The finite elements auto- 
matically provide linear inertial forces and structural 
damping, based upon the initial state, as well as nonlinear 
memees Cue to geometry and material effects. 

We will develop the basic equations for the transient 
solution. The methods were chosen to be compatible with 
existing MSC/NASTRAN static nonlinear solutions as well as 
with the existing linear transient methods. Thus, the 
Newmark Beta method for transient integration is combined 
with the modifications to Newton's method for nonlinear 
solutions. The additional iteration steps provide equilib- 
rium solutions at each time step, thereby guaranteeing 
stability and accuracy for arbitrary time step size. The 
basic equations are developed below, followed by a discus- 
Sion of the stability limits. 

a eeesLlc Equation 

Adding dynamic loads to the basic equation given 
previously, we obtain a load equilibrium error vector {6} 


at time step n by equation: 
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Cae a = Bul Ss 
where {P_}= “Average” load over the time period 
eee en Sa’ 
{u_} {U_}= Corresponding acceleration and velocity vectors. 
= "Average" elasto-plastic element total force 
VEGEOL . 
The above equation is solved at the "reduced" 
(uy) displacement vector size. The approximation errors due 
to dynamic reduction methods are not included in the error 


vector {6}. Applying Newmark's averaging method over a 


finite time period, t —ate ae the "static" forces are 
n=1 n+l 


eee} = {BF(u_ 4) + (1-25) F(u,) + 8F(u__4)} (28) 
where § = Newmark Beta operator and F(u,) 1s the nonlinear 
force due to a generalized displacement vector {U_}. An 


identical definition occurs for ay from the applied loads 
at each time step. 
From central finite differences, the acceleration 


and velocity vectors are: 


(u_} = <3 {Au_,, - du,} = — {d_} (29) 
{a} = = (Au, + du} (30) 
n Zit ntl n 
where {Au_4y = {ua - u} 
for small angles and At = thay = th is the "time step size". 
For large angle changes we require that fu}, {u}, and {Au} 
are vectors in the global coordinate direction. The vector {u} 
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contains gimbal angles requiring a tranformation R(u_) 
Such that 


{u } = tu} + ee {Au_ } (31) 


n+L1 +] 


the load error in terms of the current vectors estimates: 
iy 


a T 1 
{6} ae 2 ree) Pe ste TAL B | ce, =e B] {du} 
ope lace) Er + PE) (32) 
and 
Ee oo Lou + ine) (33) 


At any time step the vectors {P,}, {AU_}, {F,} , and {F__,} 


} and f€£ = 


are known. In general the vectors (dt, { | 


Unt 


F(u.,,) must be found either by approximating {F TL OueDy, 


n+l 
an iterative search. 

Note that for a linear solution ee equals 
laa | , Micromass themreduced stittness matrix. In this 


j 


case Eqns (32) and (33) may be used to solve for {u4y 
directly by setting {6.} = 0. With §=1/3 this will produce 
the linear SOL 27 MSC/NASTRAN results. 
Nonlinear Iterations 

For a nonlinear solution, Equns (32) and (33) may be 
solved with Newton's method (or a modified version). Using 


{d_} as the primary solution variable provides the following 


iteration algorithm. 
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First assume a linear approximation: 











{6***} ices at {a_*** - a7} = 0 (34) 
L—t 
(A, } = |2°n Peis 
n 
Then, the estimated displacement change is 
ee (ag) is cat) } (36) 


The algorithm is identical to the static case except that 
instead of a tangent stiffness matrix, the "left hand side" 


matrix, obtained from Eqns (32) and (34) is 


| il ~ 
io 
where [|K] = = is the current tangent stiffness matrix. 


At a new time step the load iterations may be started by 


1°] 
assuming that Gs rahe = {Au_t and therefore: 
(38) 


{6} = {0} 


However, to be consistent with displacements, we must extra- 


polate the nonlinear forces: 


O load = = 
tee eee eee ~ Ft (39) 


LE ey 7 


Substituting Eqns (38) and (39) into Eqn (32) for i=0 results 


maecne first estimate: 
O1_ ———_ = 
(ice ie a a [A,] {Au} ea (40) 


By 








where 


— (on (41) 


Here A is a precalculated matrix. {AU} and {FF} are 


previous results. For subsequent iterations, (1>0) Eqn (32) 
becomes: 
a) a _ ee 
Uae [Az] ido BLE Fuel? vd 
where l 1 


At each time step the code will iterate on Eqns 
(36) and (42) until (57 passes the convergence tests or 
the number of passes reaches an iteration limit. With a 
Single step, i=l, calculating only O5, the results will be 
identical to the existing NOLIN results in MSC/NASTRAN. For 
faster convergence the iterations may continue, the matrices 
ieaeeee Updated, and/or the time step size may be reduced. 
eee actors Determining Response to Transient Loading 
The response of a metal to dynamic loading depends 
upon three elements which are characteristic of the material: 
Ene dynamic constitutive relation, the pressure-temperature 
phase diagram, and the dynamic fracture criteria. In the 
absence of phase changes or fracture, the material response 
1s governed completely by the constitutive relation, or 
dynamic equation of state. The constitutive relation is 
characterized as "elastic", "viscoelastic", "visco-elastic- 


plastic", etc. It depends on macroscopic parameters such as 
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Pung s modulus and the viscosity coefficient which, ona 
microscopic scale, arise from atomic interactions and dislo- 
cation processes. The phase diagram determines whether a 
given loading history will bring the metal into a pressure- 
temperature regime where phase transitions are possible. If 
this occurs during shock loading, multiple shock waves may 
result. Microcracks and/or voids begin to appear in a metal 
under tension whenever the criteria for nucleation and growth 
of damage are met. These criteria are functions of material 
parameters and the stress history. 

In this chapter the dynamic constitutive relations 
are discussed. It is the most important element to determine 
the response of a metal to dynamic loadings. 

In macroscopic or rheological language, iron can be 
described as a viscoelastic-plastic material, l.e., it 
responds elastically until yielding (with upper and lower 
yield points), and is strongly rate-dependent. 

The constitutive relation, or dynamic equation of 
state, thus gives the stress as a function of strain and 
time, or, alternatively, can be written as a differential 
equation in OeGre and an important goal of present day 
research is to develop the capability to derive this macro- 
scopic constitutive relation by averaging the effects of 
microscopic processes. 

This is analogous to the way in which the thermo- 


dynamic equation of state of a gas is derived by averaging 
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the effects of the molecular collisions. Once developed, 
this capability would be an extremely powerful tool, since 
the response of a metal to dynamic loads could be predicted 
from Known microscopic material properties. Although this 
goal has not been attained as yet, there has been significant 
progress during the last few years in this direction for the 
case of crystalline materials, in which dislocation processes 
play the role of the “molecular collisions". In iron, the 
density of mobile dislocations is of the order of Lo, ene, 
so the number of dislocations present is certainly large 
enough to make a statistical approach valid. 


The basic approach is based on the relation from 


dislocation theory: [Ref. 3] 


E 


Ve NO bv (44) 
mere Yy is the plastic shear strain rate 
"ae ee | | | 
cee = 2 ‘Sx, ~ ox, ) , where ¢, is the displacement in 


the x. direction), No is the mobile disclocation density, 

b is the Burgers vector, and v is the average disclocation 
velocity. This equation provides the link between microscopic 
dislocation motion and macroscopic plastic straint rate, and 
is based on the fact that it has been found possible for many 
crystalline solids to express NO Sows yas functions of y* 


and applied macrosopic shear stress, T. 
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The constitutive relation is obtained by combining 


Eqn (44) with the incremental Hooke's law 


e e 
= +dAd — 
one oe ae ExK one (25) 


under the assumption that the total strain is the sum of 


elastic and plastic components: 


€ = ¢~ eF 
1) iy 3 hg (46) 
where XA and nh are the Lame constants, 6.. is the 


1) 
| Kronecker delta, and a repeated index indicates summation. 





The plastic strain is assumed to cause no volume change, i.e., 
cf = 0 (47) 


The combination of Eqns (46)-(48) yields, after taking the 
time derivative: 
212 


ee -2 Weis 7 ere ae = 2 are 


(43) 
In order to combine Eqn (48) with Eqn (44), we must 
First specify the loading geometry. In plate-slap tests 


under conditions of uniaxial strain, all total strain com- 


ponents except 25 in the direction of wave propagation are 
zero, and Ex Ey: Then equation (48) becomes 
pa (he2uven = -2 wet (49) 
x ok x 


Also, in this case, the plane of maximum shear stress in 
isotropic media is inclined at 45° to the direction of shock 
Beepagation, and 


Joe PS. 3 IP 
Yp =7 f€, - &€ y) FTE (50) 











Tm dynamic constitutive relation for this loading geometry 


mes OOtCained by combining Eqn (44), (49), and (50): 


: a . koi 
oy (A - 2u)e, = 5 HN bv GSar) 


iMeeeeA, ww, and the right-hand side of Eqn (51) are known 
as functions of macroscopic stress and strain, then the 
equation can, in principle, be solved for any problem 
Bavoelving Uniaxial strain. 
 Lynamic Plastic Deformation of Plates 

We consider cases of the impulsive loading of thin 
metal plates, a simple situation in which an exponential 
shock wave impinges on a clamped square plate. In this case, 
the plastic deformation is supposed to be brought about by 
bending only and combined bending and tension is not analy- 
tically considered. The pressure in the shock wave may be 
assumed to be given by [Ref. 4]. 

P =p, exp(-t/t) 

where p is the peak pressure at a point at the head of the 
wave, t denotes time and tT is a time constant --- the time 
taken for the pressure to fall to aa Pislael al “seepage g\atle ial 
explosives is ~30 secs. We obtain simply, some of the 
results of Cox and Morland [Ref. 3], who theoretically 
investigated the dynamic deformation of plates, side length 
2a, position-fixed at their periphery, when subject to a 
uniform constant transverse pressure p for time Tt. Their 


detailed examination showed that two categories of load may 
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be distinguished; medium load when Pas a) eS 2p. and high 
load when 2p. p, where Pe denotes the pressure required 
mmeeeto Gnitorce plastic collapse of the plate, sittatically 


For p we have, see Fig. 3-l. 


dav + W/2- ML, = 5: 4a Sap oe (52) 


- 2 
1.e. P = oM,/a 


For simplicity we deal only with their first case, 
and for the second, more complicated one, the reader 1s 
referred to the original source. For this medium range of 
loading two sub-divisions are required: (1) when 0<t<tT 
and (2) when t>T. Use is made of the results obtained 
above for static situations using hinge lines, and we take 
over whatever results may be useful. As previously, the 
analysis pays no special attention to material rate effects. 
Purther, it 1s implicit in this approach that the static 


and dynamic modes of deformation are identical. 


FRONT NET ON THE 


LOAD 
LOAD 





Uy TIME 
TIME 


Prieute 3,1. Wynamic Loading of Square Plate. 
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a. Medium Load: O0O.< t < T 

Assume the mode of deformation of the plate is 
that of Fig. 3.2 and that the plate diagonals become hinge 
lines. Since p is constant and greater than Pos and 
because the four equal parts of the plate act as rigid bodies 
each rotating about its boundary or side, they must undergo 
an angular acceleration, a. If changes in plate geometry 
throughout the deformation are neglected, the resistance to 
deformation, which arises from the yielding taking place in 
the hinge lines, will be constant; this resistance is related 
to Pos the pressure to cause static yield. Thus the excess 
load is constant and hence a must be constant. The work 
done by the excess load will manifest itself as kinetic 
energy acquired by the four rigied regions of the plate, l.e. 
4 x = Mmiyeewhere = is the moment of inertia of a triangle 
about a side of the plate and w is its current angular 
velocity. Because a is constant, at the end of time the 
plate centre will have descended through a distance vt/2, 
where v 1s i1tS current speed. Thus, an equation involving 
the work available for giving rise to rotational kinetic 


omemgy OL the plate triangles is, 





ey 
(p> - p_) 4a*° A eee (53) 
S 3 2 
Hence substituting for Pe PON Mee enOtingsthat I = ma*/6 
ama Simplifying, 
6M (p-p_)t 
(p - —B yt = 772 or v = 2 ——=— (54) 
a 








or directly by momentum consideration, thus, 





4 
., . a2aalma ly 
((p Pp) t 5 ) 5 7 = (25) 
The deflection Wy ae the end of time TT is, 
a /2 6 (p - p_)t*/m (56) 


(2) Medium Load: T>t 

After the removal of pressure p at t= tT, further deflec-~ 

Eon, Wor is added to Wie because the rotational kinetic 

energy of the four plate segments must be dissipated in doing 
plastic work in the hinges. Let the segments rotate through 


further angle o> before coming to rest, then, 





2 W 
4 - (ee ee 4 av¥2 - M ($72) = 8aM - Af 5 
2 a7 D p 
Hence 
22 
= ma*v*  _ my“ = ee (57) 
Z 2 4 M, 4p. Mm De 


using Eqn (55). Thus the total central deflection is, 


7 p,) t°/mp (58) 


= 
7 
= 
t- 
+ 
= 
t 
Te 
oO 


If the total response time of the plate is note that 


The 
T,P. = pt. In the case of thin plates clamped at their outer 
periphery it is clear that after a significant amount of 
deflection there must be stretching. An analysis of the 
dishing of peripherally clamped metal plates when given a 


transverse initial speed of a few hundred feet per second 


has been given by Hudson [Ref. 3]. Descriptively, Hudson's 


a5 








model is easily understood by reference to Fig. 3-2. Central 
portions of the plate maintain their given transverse speed 
until overtaken by a bending wave propagated radially inwards 
from the clamping ring. Fig. 3-2(a) shows the imagined plate 
configuration at successive intervals during the deformation 
process; when the plate deflection is significant, the process 
1s predominantly one of stretching and for many practical 
purposes the technological approach is adequate. 

Figure 3-2(a) shows the terminal state taken up 
by a thin rectangular steel plate clamped around its periphery 
after being subjected to the uniform but intense impulsive 
loading, the maximum speed of the plate is a few hundred 
feet per second, but contracting in area with time; the 
dynamic situation during deformation, with plastic hinges 
moving into the plate from the clamped edges and intersecting 
along lines equally inclined to the sides from the corners, 
1s indicated diagrammatically in Fig. 3-2(b). 

Now the yield criteria of metals should be con- 
Sidered. Two approaches to the interpretation of a yield 
criterion are possible; one is Tresca's criterion, the other 
is Von Mises criterion. But for most metals it 1s generally 
found that while the experimental points fall between the two 
ellipses, they incline towards Von Mises' ellipse. Experiments 
were earlier performed by Lode (1925) on tubes subject to 
tension and internal pressure. Bending and torsion was used 
by Siebel in 1953. All results confirm the general opinion 
that the best simple yield criterion for metals is that of 


Von Mises [Ref. 5]. 46 
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Figure 3.2. Dynamic Deformation of Fixed Square Plate. 
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Be PROBLEM DESCRIPTION AND NASTRAN DECK SETUP 
mm Peepnem Description 
The mode is shown in Fig. 3.3. It is a 20 in. by 
20 in. steel square plate model which is divided to 25 QUAD4 
(NASTRAN card) elements and its boundary is clamped. Its 


thickness is 0.8 inches and the engineering data as follows. 


Ee— > .0Bt? 

v = 0.3 

ma §20..0- 1m. 

9 = 7.33E-4 lb/in? 


moe O20 in. 
where E is the elastic modulus, v is the Poisson's ratio, 
Lis the length of the plate, p is the mass density, and 
t is a thickness of the plate. 

Next the shock wave applied impinges on the steel 
plate in the -Z direction at time 0.0 seconds (Fig. 3.4) and 
its peak overpressure is 3.0E+4 psi and the shock wave dura- 
tion time (ty) is 5.0E-3 secs. To find the pressure, the 
following formula is used [Ref. 6] 

p= p,(l-t/t,) exp(-at/t,) 
where Pan 1s peak reflected overpressure at time 0.0 secs 
and a 1S wave form parameter. But the shock wave pressure 
cannot be applied as a load to the plate since the front face 
shock pressure of the plate is relieved by the rear face shock 
wave which is positive Z direction. The following formula is 


used [Ref.4]}. 
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p = P exp (-t/T) 


The load data as follows. 
peak lateral loads (-2Z direction) 
So = 3.Bt4e9si on the Elenent (213) 


2.2E+4 psi on the Element (207, 208, 209, 2ie2, 


P 
Ze 2 ao, 219) 


Von Mises criterion would be applied for this model. The 


time step (t) to be applied is 1.0E-6 secs from 0.0 sec to 


7 .0E-4 secs. 
/ 


20500 20501 yAOho) 0 20503 20504 20505 





(201) (202) 


20QQ0 IOGL OOO 


Figure 3.3. Modeled Steel Square Plate. 
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Figure 3.4 Shock Wave Applied. 
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Figure 3.5. Loading History. 
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2. NASTRAN Input Data Deck 


Reference is made to the echo of the input data deck 


which appears in the output. 


10 
a 
ibe 
i 
14 
igs 
16 
ay 
18 
ing 
20 
ii 
22 


Zz 


a. Sol 66 [Ref. 6] 
NASTRAN PREFOPT=2 
ID NONLINEAR, STATIC 
SOL 66 
TIME 30 
CEND 
TITLE-NONLINEAR DEFORMATION WITH STATIC LOADS 
SUBTITLE=TRANSIENT RESPONSE BY SHOCK WAVE LOAD. 
ECHO=UNSORT 
DIS P=ALL 
QmOAD=ALL 
SPCF=ALL 
Ger OR=ALL 
ELST=ALL 
See L-ALL 
SPC=200 
SUBCASE 1 
HOAD=10 
NLPARM=10 
BEGIN BULK 
$ 
S$ LATERALLY LOADED FLAT PLATE 


S 


oul 








24 
25 
26 
ot 
28 
29 
30 
Bt 
32 
33 
34 
3) 5) 
36 
oy 
38 
2 
40 
41 
42 
43 
44 
45 
46 
47 
48 
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MESHOPT,,,,,YES 

$ 

S$ USE MSGMESH TO GENERATE 5*5 ARRAY OF QUAD4 ELEMENTS 
> 

EGRID,1,,0.,0. 

MentD,2,,20.,0. 

HGRID,3,,20.,20. 

BERTD,4,,0.,20. 

OpTDG,2,,,5;71,-2,-3;,,tFuD2 

+FLD2,5,-4 

CGEN,QUAD4,201,200,2 

S$ REFER TO PAGE 2.18-7,APPLICATION MANUAL. 
PSHELL,200,65,0.8, 65,,65 
MAT1,65,3.0E+7,,.3,7.33E-4 
MATS1,65,40,PLASTIC,,1,1,2.6E+4 

TABLES 1,40,0,,,,,,,+TABL 
+TAB1,0.0,0.0,0.001,3.0E+4,0.025,4.-E+4,0.05,5.8E+4,+TAB2 
+TAB2,0.10,6.7E+4,0.125,7.0E+4,0.2,7.5E+4,ENDT 
PARAM, COUPMASS ,1 

PARAM, K6ROT,1.0 

PARAM,LGDISP,1l 
PeenD2,10,=-14.7,208,212,213,214,218 
NLPARM,10,1, ,AUTO,,20,W,YES 
SNLPARM,20,3,,AUTO, ,20,W,YES 

SPCG, 200 ,2,123456,A,B 


peeG,200,2,125456,B8,C 
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20 
SL 


DZ 


opeG,200,2,123456,C,D 
Seee,200 ,2,123456,A,D 
ENDDATA 
1) Executive control card 
card 1 - The NASTRAN card contains the keyword PREFOPT 
= 2 which requests that MSGMESH be used to 
generate input data. 
card 3 - Request nonlinear static analysis, Rigid format 
66. 
@e Case control card 
card 15 - Requests superelement generation and assembly 
* should be set as ‘all' 
card 18 - Requests static load combination. 
* SID should be matched to the PLOAD2 card in 
the bulk data deck. 
card 19 - Requests nonlinear parameter selection. 
* SID should be matched to the NLPARM card 
in the bulk data deck. 
* This card may appear above or within a 
subcase. 
* LOADcard is only applicable in statics, 
buckling, heat transfer problems. 
3) Bulk data cards 
card 24 - Requests out options control 
* An entry of YES causes the printing of the 


Prmoiemecatad echo  . 


iS 





card 28-3F - Defines vertices of quadrilateral field 


Cawed 


card 


card 


card 


card 


card 


card 


card 


55 


34 


38 


Sy, 


43 


44 


45 


46 


using MSGMESH EGRID command. 
Defines a field of grid points and generate 
Grave Candas . 
* The grid number generated by MSGMESH starts 
fimemethe FID Of GREDG (if FID is 2, the first 
@ria Doint ID is 20000) 


* 5 means 5 elements in each direction. 


- Generate element connection. 


* QUAD4 is the element type which should be 
generated. 

* Element ID number starts from 201. 

* Fifth field 200 is the PID of pshell card. 
Specifies table references and material 
properties which are stress-dependent for use 
in material nonlinearity applications 

* Taken Von Mises yield function, isotropic 
hardening rule. 

Represent stress-strain data. 

The value 1.0 suppresses singularities 

* It is intended primarily for geometric 
nonlinear analysis. 

The value 1 allows the nonlinear calculations. 
Request the initial condition on the elements 
in the negative Z direction. 


Controls the nonlinear analysis iteration. 
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10 
tal 
12 
ES 
14 
LS 
16 
ey 
18 


Eo 


x To re@tiee the ouleput, use sixth field. 
@ er Gapid CoNmvergence, increase the Error 
tolerance and maximum iteration and the limit 
on diverging iterations. 
card 48-5l - Generate spcl cards. 
This card is only used to MSGMESH generation. 
bee SOl 99 
NASTRAN PREFOPT=2 
ID NONLINEAR, TRANSIENT 
Som 99 
Tae 59 
CEND 
TITLE=NONLINEAR DEFORMATION WITH DYNAMIC LOADS 
SUBTITLE=TRANSIENT RESPONSE BY SHOCK WAVE LOAD. 
ECHO=SORT 
Ser 1=20202, 20203,20204,20304,20404 
ee — AO), 206,212,213,214,219,225 
few 5=20000,20003,20005, 20200, 20400 
SDISP=1 
OLOAD=1 
SPEF=3 
EES T=2 
SEALL=ALL 
SPC=200 
LOADSET=77 


SUBCASE l 
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20 
aoks 
I 
23 
24 


Zo 


26. 


27 
28 
29 
30 
Sidi 
32 
338 
34 
5) 
36 
37 
38 
ae 
40 
41 
42 
43 
44 


45 


DLOAD=37 

TSTEPNL=45 

SUBCASE 2 

DEOAD—37 

TSTEPNL=46 

SUBCASE 3 

BEOAD=s 7 

TSTEPNL=47 

OUTPUT ( XYPLOT) 

PLOTTER NAST 

CSCALE 1.3 

XTITLE=TIME IN SECS 

XGRID LINES=YES 

YGRID LINES=YES 
YTITLE=DISPLACEMENT GRID 20202 
XYPLOT SDISP RESP/20202(T3) 
Ween E—VELOCITY GRID 20202 
XYPLOT SVELO RESP/20202(T3) 
YTITLE=DISPLACEMENT GRID 20101 20102 20202 
XYPLOT SDISP RESP/20101(T3) ,20102(T3) ,20202(T3) 
BEGIN BULK 

$ 

S$ LATERALLY LOADED FLAT PLATE 
$ 

MESHOPT,,,,,YES 


° 
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46 
47 
48 
49 
10) 
ou 
a 
55 
54 
oe 
56 
mj 
58 
5, 
60 
ol 
62 
6 3 
64 
65 
66 
67 
68 
69 
70 


7k 


S USE MSGMESH TO GENERATE 5*5 ARRAY OF QUAD4 ELEMENTS 
$ 

BGrip,l,,0.,0. 

EGRTD,2,,20.,0. 

MeriD,3,,20.,20. 

BERTD,4,,0.,20. 

eG, 2, ,)5,-L,)-2,-3,,+FLD2 

+FLD2,5,-4 

CGEN,QUAD4,201,200,2 

S REFER TO PAGE 2.18-7,APPLICATION MANUAL 
PSHELL, 200,65,0.8,65,,65 

Merl, 65,3-0B+7,,.3,7-338-4,,,0.01 
Mersl,65,40,PLASTIC,,1,1,2.6E+4 
TABLES1,40,0,,,,77,,+TABL 
+TAB1L,0.0,0.0,0.001,3.0E+4,0.025,4.0E+4,0.05,5.8E+4,+TAB2 
+TAB2,0.10,6.7E+4,0.125,7.0E+4,0.2,7.5E+4,ENDT 
PARAM, COUPMASS , 1 

SPARAM, AUTOSPC, YES 

PARAM,W4,1.616E+3 

PARAM, K6ROT,1.0 

PARAM, LGDISP,1l 

PARAM, SLOOPID,1 

PARAM, LOOPID,1 

PARAM, STIME,0.000450 

PARAM ,MAXLP, 7 


PEON oe, Lure. 0,7,L.0,6,1.0,5 


=| 








Te 
oS 
74 
75 
76 
77 
78 
re) 
80 
eb 
82 
83 
84 
85 
86 
87 
88 
89 
210 
OE 
OZ 
93 
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Ord, 7, Gden,.0.,.22 

mam®,77,61,71 

MLOAD2, 71, 3.2E#4,213 

REOADL, 6,62, ,0,22 

omg, 77,62,81 

BmeAD, Sl,2.56+4,207,208,209,212,214 
mEOAD!,5,63,,0;,22 

momo, 77,63,91 

mon 2,91,2.5E+4,217,218,219 
TABLEDL,22,,,-777,tNEXTL 

foe, ~l.0-0,0-.0,-0.0010,0.0001,-1.0,0.0010,-0.560,+NEXT2 
+NEXT2Z,0.004,-0.001,ENDT 
RereePNL,45,150,1.0E-6,20,AUTO,,20,W,+LL1 
fm, 1.0E-7,1.0E=-7,1.0E-7,20,20,10 
TSTEPNL,46,200,1.0E-6,20,AUTO, ,20,W,+LL2 
fame, es OER 5, 1.08-5,1.0E-5,20,20,10 

Rome PNL,47;,120,3.0E-6,20,AUTO,,20,W,+LL2 
goa, 1.0E—4,1.0E-4, 1.0E-4,20,20,10 
eeeG,200,2,123456,A,B 

meG, 200,2,123456,B8,C 

omeG, 200,2,123456,C,D 
SPCG,200,2,123456,A,D 

ENDDATA 

1) Executive control cards 

card 3 : Requests nonlinear transient analysis, Rigid 


format 99. 


Sys. 





2) Case control cards 

card 18: Selects a static load set for use in dynamics. 
Thus this may be referenced by dynamic loads. 

card 20: This card selects a dynamic load for nonlinear 
transient problem. 

card 22-24: Requests restart at time 0.00015 secs. 

card 25-27: Requests restart at time 0.00035 secs. 

card 28-39: Requests the xy plotting which is 

corresponding to the time step. 

3) Bulk data cards 

card 64: This gives damping effects for transient 
analysis. 

* Sometimes, more rapid convergence effects 
can be obtained by this damping factor with 
damping coefficient in MAT1 card (GE). 

card 67: This card identifies the initial conditions 
from a previous SOL 66 nonlinear static 

So LuEve@n . 

* Refer to the param card, SLOOPID. 

card 68,69: Request these two cards when we taken more 
than one subcase. 
card 70: Requests more internal calculation. 
* We can reduce iteration steps and prevents 


the overprinting due to iteration steps. 


5) 








Gama Jl: Defines: a dynamac loading condition for 
transient response problems as a linear 
combination of load sets 
* See the remarks 5, 7, 9 of DLOAD card in 

bulk data. 

card 72, 75, 78: Defines time dependent dynamic load. 
* Requests a static load card set (LESQ). 

eommec 73, /6, 1/9: Defines a penGence of static load sets. 

This load sets may be referenced by dynamic 
load card (TLOAD1) 

* This card should be used when we want to 
apply area pressure on elements than a point 
pressure on grid point. 

card 74, 77, 80: Defines a uniform static pressure load 

on the plate elements. 

card 81-83: Defines a tabular function for use in 

generating a time-dependent dynamic loads. 

* In this card the peak pressure (3.0E+4 psi) 
occurred at 0.0001 sec. 

card 84, 86, 88: This card control the nonlinear 

transient analysis. 

* This card is expalined in the sensitivity 


analysis. 


eee oENSITIVITY ANALYSIS 
1. General 
In transient analysis, two types of instability could 
Scelmethewtirst is the familiar nonlinear load ateration 


divergence which also occurs in static analysis. The second 
60 





is the divergence which grows with time in the transient 
integration. Both instabilities are caused by uncorrected 
nonlinear equilibrium errors. 

A convenient method, from Von Neumann, for analyzing 
the stability limits is to assume that the nonlinear forces 
are nearly linear and the error vector has a constant conver- 
gence. It is assumed that the nonlinearity has a first order 


approximation: [Ref. 7] 


1a N 
{F(u)} ~ (KX + ax™?] (u} (59) 
where axNt is the difference between the tangent stiffness 
matrix and its approximation kK’. The error vectors, {6}, 


are assumed to grow at the rate 4A, defined as: 
{s***}= 016.) (60) 


Note that if |A| > 1 the system will be defined as unstable. 
After lengthy calculations we may summarize the various 
Criteria, assuming that the matrices are reduced to equivalent 
scalar modal quantities, the criterian for stable solutions 
are: 
1) For time steps with converged "static" iterations: 
ieee 6 2 1/4 
This is the same criteria as linear analysis. 


2) For "static load" iterations: 


B A eo Sas KT + —_ ee — 
At 


Dat _ 
This states that the mass and damping add to the effective 


linear stiffness fe? | and improve the stability. 


ou: 





3) For time step integration with no intermediate static 


iterations, as in the standard Newmark Beta method: 


Wee = (46 -1) x + 


M 
At? 





This restriction 1s more severe than criteria b above, 
proving that the internal iterations are more stable 
than the Newmark integration. 

In summary, the Von Neumann method will have fewer 
divergence problems than either the static nonlinear solu- 
tion or the single step transient nonlinear methods. The 
better stability of the method and the capability to use 
larger time steps outweighs the cost of a few internal 
iterations on the static element forces. 

2. Nonlinear Transient Controls 

The input data required for SOL 99 is a combination 
of transient control data, similar to SOL 69 (Direct Transient, 
superelements), and nonlinear modeling data similar to SOL 66 
(Nonlinear Statics). The nonlinear properties are defined by 
nonlinear material data (MATS1), gap elements (GAP), and 
large displacement effects (parameter LGDISP). The transient 
effects are produced by loading functions (TLOADi, DAREA, etc.), 
damping (parameters, elements, and material data), and element 
mass properties. 

The only unique data required for SOL 99 is supplied 
on the TSTEPNL data cards. (Several optional parameters are 
provided for tuning the method and for use in restarts.) The 
TSTEPNL card in itself is an amalgam of the normal transient 


TSTEP card and the nonlinear control card, NLPARM. 
ay 





a. Cage Control 
che gughatiicane, cules for nonlamear transient are: 
1) Each subcase defines a time interval starting from the 
last time step of the previous subcase which is subdivided 
meecO smell time steps. The output TIME printout will be 
labeled by the net time, including all previous subcases. 
2) The input loading functions may be changed for each sub- 
case or continued by repeating the same DLOAD request. 
The bulk functions are evaluated at the net model time, 


Which starts at the last time step of the previous subcase. 





3) Each subcase may have a different time step size, time 
interval, a Piemataon Control selected by the TSTEPNL 
request. The case control requests which may not be 
Ghanged after the first subcase are: SPC, MPC, DMIG, and 
ae. 

4) Output requests for each subcase are processed indepen- 
dently and the output for each case are generated before 
the next subcase is processed. Separate XYPLOT ouputs will 
occur for each subcase. Available outputs are DISPLACEMENT, 
Pemsocrhy, ACCELERATION, OLOAD, STRESS, FORCE, SDISPLACEMENT, 
SVELOCITY, and, SACCELERATION. Note that NONLINEAR and 
SPCFORCE outputs are not provided. 

Aimebttial conditions are not explicitly provided for user 
input in SOL 99. The reason is that the nonlinear element 
configuration and equilibrium state must be known at every 
time step. If initial conditions were given, all of the 
nonlinear element forces and plastic stresses, as well as 
the past history, must be given. Therefore, the initial 
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conditions become part of the analysis and are generated 
by initial transient subcase or by restarts from a previous 
Seoelc solution. 

6) All normal superelement model generation options and 
matrix reduction options are allowed for the linear portion 
of the structure. General Dynamic Reduction, Component 
Mode Synthesis, and Guyan Reduction may be performed for 
nonresidual superelements. Therefore, the residual super- 
element may contain scalar degrees of freedom representing 
linear model formulations. 

b. TSTEPNL Data 

The TSTEPNL data card description is provided in 
Section 2 of the MSC/NASTRAN User's Manual. The input fields 
Specify the time step size, the number of steps, and the out- 
put interval as well as the nonlinear iteration options. 

The nonlinear iteration data are provided for 
controlling the accuracy and stability of the nonlinear 
calculations and they will not effect linear solutions. The 
controls are similar to those on the NLPARM bulk data used in 
Static nonlinear analysis. 

As discussed above two types of errors or insta- 
bility are encountered in the transient solution. The load 
iterations are similar to static load analysis and may diverge 
at one or more time steps when large nonlinear stiffness 
changes occur. Load iteration divergence may be cured by 


smaller time steps or by stiffness matrix updates. The 
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diagnostic output value LAMBDA-I (from DIAG 50) indicates 
mee Tate of growth of this error. 

The second type of divergence occurs in the time 
domain when the load iterations do not converge, or are limi- 
ted by the user. The errors may increase for each new time 
step until they cause the code to abort. This problem may be 
cured by allowing more load iterations or by using smaller 
time steps. The diagnostic output value LAMBDA-T (from DIAG 
50) indicates the rate of growth of this error. 

The size of the time step is of primary importance 
in transient analysis to control accuracy and avoid diverging 
solutions. The use of several subcases with different time 
step sizes to correspond with nonlinear activity is recommended 
as a normal procedure. When the user wishes to avoid overly- 
small time steps for highly nonlinear zones, he may increase 
the iteration limits and/or the error tolerance and the diver- 
gence limits. 

c. Restarts and Data Base 

Solutions 66 and 99 share the same data base/ 
Superelement file storage logic for the nonlinear tables and 
matrices. Therefore, the restart system for transient analysis 
is allowed to use either a previous static or transient non- 
linear analysis as initial conditions. 

For restarting from previous static analysis only 
the first subcase is affected. Simply provide a data base 


created in SOL 66 and specify the parameter. 
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PARAM,SLOPPID,N 
where N was the printed value of LOOP for the desired 
S@atic solution. Constraint sets, direct input matrices, 
mass, and damping may be changed. 

Restarts from a previous transient execution are 
available for a number of cases. If the same model is to be 
reexecuted, only the residual superelement needs to be re- 
assembled (SEMA, SELA=0). If the end results for the previous 
transient run are to be used for the initial conditions at 
t=0, add the parameter LOOPID=N where N corresponds to the 
Nth subcase of the previous run and N dummy SUBCASE commands 
to start the residual case control execution. 

The normal restart for a transient run is to 
continue from the last step of a previous subcase with 
different loads and/or TSTEPNL data. Provide the following 
Parameters: 

LOOPID = N - Start from the Nth subcase 
Stime = t - Start from time t. 
Note that constraint sets should not be changed 


to avoid incompatible matrix sizes. 


Dene NUMERICAL RESULTS 

The progressive deformations of the plate are plotted 
along the x-axis at five different time locations as shown 
Mmmeerog. 3.6(65). At the early tyme (1.5E-4 secs), the response 
of the plate is in the elastic region. During the deformation, 


plastic hinge lines are formed along the clamped edges and the 


66 








G@iageonals of the plate as indicated in Fig. 3.6(a). This 
response phenomenon agrees well with Hudson's studies 
[Ref. 4]. As described in Section A-4 of Chapter 3, when 
the plate deflection is significant, the plate experiences 
in-plane stretching. 

Time history displacement responses of Grid points 
Peeeees, 20201, and 20202 are plotted in Fig. 3.7. As shown 
in the figure, for the short duration impulse, the maximum 
displacement occurs at the center of the plate after termi- 
nation of shock loading. The rate of increase in displacement 
is the maximum at the termination of shock loading as shown 
mameero. 5.5. During the application of the shock loading, 
the acceleration level is high, implying the high inertia 
force in the system. At the termination of the inertia 
force in the system, the acceleration stays in the constant 
level until the plate displacement reaches the maximum. 

Time hisotry Von-Mises stress of QUAD4 element 213 (plate 
center) is plotted in Fig. 3.7. The stress variation is 


highly non-linear due to the material non-linearity. 
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Figure 3.6(a). Plastic Hinge Line Formation. 
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(b). Progressive Deformation of the Plate 
Along Section A-A. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The state-of-the-art literature review revealed that the 
skirt plate is very effective to absorb initial explosion 
energy before the shock loading reaches the main armor. The 
material of the armor must possess the high hardness to prevent 
the severe damage and the high ductility to minimize the 
spalling. The composite materials may serve more effectively 
than the conventional material to minimize the damages against 
the powerful modern projectiles. 

The non-linear finite element analysis of the plate sub- 
jected to the air shock loading verifies the plastic hinge 
line formation and in-plane stretching inducing snap-through 
collapse. For the short duration shock loading, the maximum 
deformation occurs after the termination of loading, and the 
velocity reaches the maximum at the termination of shock. 

The Von-Mises stress is highly non-linear due to the material 
non-linearity. For the non-linear analysis using MSC/NASTRAN, 
TSTEPNL bulk data card is found to be the most important card 
which controls the accuracy and the stability of non-linear 
calculation. 

To insure the convergence of solutions, convergence limit 
in the TSTEPNL card in the bulk data deck should be considered. 
If a program was run with big convergence limit, the solution 


Will be gotten rapidly with small iteration steps but it is 


The 





not accurate and large time steps can't be taken. On the 
contrary, if the program was run with small convergence 


limit, the solution will be accurate but it is gotten slowly 


with large iteration steps. Thus, it is expensive and anyway 
large time steps can't be taken too. Thus, compromise should 
be taken to run a large time steps. For example, if one non- 


linear transient problem is needed total 700 time steps to get 
the solutions, the following procedure is recommended. In the 
first subcase, the program is iterated with smaller conver- 
gence limit, 10*E-7 to get more accurate solution and in the 
second subcase, the program is run with limit, l10O*E-5, and in 
the third subcase, the program is iterated with convergence 
ieomese, lO*E-4. 

Consequently, the total 700 time step will be calculated 
Successively. 

The sensitivity studies on shock period and amplitude is 


suggested to understand their effects on plate responses. 
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APRPENDIA A 
THE COMMAND TO OPERATE DATA BASE 
owetne COMMAND TO ERASE THE DATA BASE NAMED "MSS .S2770P .NODE" 


//ERASE JOB (2770,0252) ,NASTRAN,CLASS=A 
//*MAIN ORG=NPGVML.2770P 

// EXEC PGM=IEFER14 

//SYSPRINT DD SYSOUT=A 

//DD1L DD DSN=MSS .S2770.NODE,DISP=(OLD, DELETE) 


Ve 


Peete COMMAND TO FLOT FROM THE DATABASE NAMED "MSS.S2/70.PPLOT" 


/ /LEEPLCT JOB (2770,0252) ,NASTRAN ,CLASS=C 
//*MAIN ORG=NPGVML.2770P 

//*FORMAT PR,DDNAME=SYSVECTR, DEST=LOCAL 

// EXEC NASTPLOT,PLTDSN='MSS.S2770.PPLOT' 
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APPENDIX B 


DATA DECK FOR PROGRAM 


GENERAL DESCRIPTION OF OATA DECK 
The input deck begins with the required resident operating system control cards. The type 
and nuuber of these cards will vary with the installation. Instructions for the preparation of 
these control cards are given in Section 7.6 of the MSC/NASTRAN Application Manual. 
The operating system control cards are followed by the MSC/NASTRAN Data Deck, which consists 


of the following three sections: 


Executive Control Oeck 
2. Case Control Oeck 


2 Bulk Data Oeck 7 


The Executive Control feck and the Case Control Oeck both have free-fieid formats. Only 
colugns 1 through 72 are used for data. Any {nformation {n columns 73 through 30 may appear ifn 
the printed echo, but the data will not be used by the program. As explained {in Section 2.4.1, 
limited use {s made of the data in columns 73 through 80 for the Bulk Oata Oeck. I[f the last 
character on a card {s a comma (not necessarily {in column 72), the next card {s a continuation of 
this physical card. Any numoer of continuation cards may be specified, and together they forma 
logical card. 

The NASTRAN card {s used to change the default values for certain operational parimeters, 
such as buffer size or the number of data lines printed per page. More than one NASTRAN card may 
be present. NASTRAN cards are optional, but, {f present, they must be the first cards of the data 
deck. The NASTRAN card {is a free-field card (similar to cards {nm the Executive Control Deck). 


Its format {s as follows: 
NASTRAN keyword, = value, keyword, » value, ... 
An alternate format {s: 
NASTRAN keyword, = value 


NASTRAN keywords = value P 
is 





Input Data Card EGRIOD Grid Point Position 


Description: 


Defines the position of a geometric grid point of the structural model 


Format and Examole: 


X1 ,X2 ,X3 


Remarks: 1. 


Contents 





Grid point identification number (Integer > 0) 


Identification number of coord{nate system in which the position of the grid point 


is defined (Integer > 0 or blank) 


Position of the grid point in coordinate system CP (Real or Blank) 


The meaning of Xl, X2, and X3 depends on the type of coordinate system, CP, as 
follows: (see C3RDi card descriptions) 


Rectangular 
Cylindrical 3(degrees) 
Spherical a (degrees) 4(degrees) 





The EGRID card is used in MSG4ESH to define curved edges and for other special 


Purposes (see Part I, Section 3.0). 
will not be processed by MSC/NASTRAN. 


It is similar to a GRID card, exceot that it 


Entries in fields 7, 8 and 9 are not checked by MSGMESH. 
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Input Data Card GRIOG Linear and H{gher Order Interpolation 
Description: Defines a field of grid points and generates GRID cards 


Format and Examole: 


10 
GRIDG CD/CI SELB 
—— “1 61 
| 
re il 
a - 
| |sencag) | ez(ag) | en(ac) | ezcscy | ecco) | ea(co) | excoa) | eztoa) 
m |] | {| | {e{s|[«|r7_ 
| F1(EF) | €2( EF) EL(FG) | £2(FG) EL(GH) | €2(GH) | EL(HE) | E2(HE) 
: EL( AE) E2( AE) EhCGr) 1) E2ter) E1( CS) E2( CG) E1(0H) E2( 0H) | 
Field Contents 
FID [D number of a GRIDG field (Integer, 1 < FID < 99). Must be unique for al! GRIDG 
cards. when entry {s negative, internal EGRIO cards are generated. 
€0/ Cl Prescribes displacement and/or interpolation coordinate systems (Integer or 
blank). See Remarks 4 and 5S. 
PS PS entry on generated GRID cards (any of the digits 1-6 with no imbedded blanks) 
(Integer > 0 or blank) 
bent st8 Grid field dimensions when an entry, L, M, and/or Nis positive. (See Ficure 3.) 
Absolute value is a LID entry on a LIST input card, when an entry is negative. 
Se50 ScID entry on senerated GRID cards (Integer > 0 or blank) 
Gay. ..,GH Grid point identification numbers, GIOs, of grid field vertex points (Integer # 0 or 
blank). See Remarks 3 and 6. 
er CP entry on generated GRID cards (Integer > 0 or blank), and grid field inter- 
polation coordinate system definition. See Remarks 4 and 5. 
ELCXY) ID number of edge grid point (closer to X) on side XY (Integer # 9 or blank). See 
Remarks 3 and 6. 
E2(XY) ID number of edge grid point (closer to Y) on side XY (Integer # 0 or blank). See 


Remarks 3 and §. 


(Continued) 


Tq 





Remarks: il. 


GAIGE (Cont. ) . 


See Part I, Sections 2 and 3 for discussions on the GRIOG fnput card. 


Required data by type of field 


Fel ies cos 








PhO eaeGA, GG, aC 
Pll Veet.) Ga, os, GC, GO 






Freee il, Gn thru GH 


Only the parent and “(first continuation card are applicable for linear 
interpolation. 


Vertex and edge points are defined on GRID or EGRID input cards, or on GRID cards 
generated by previously processed GRIDG f{nput cards. Field vertex points are 
automatically equivalenced to corresponding entries GA, ..., GH, when these entries 
are positive. GRID cards are generated only for active field grid points. See 
Part -l,. seccron 0.1. 


CO/CI entry may be used to prescribe both the displacement coordinate system, CO, 
and the fnterpolation coordinate system, CI. Both CD and CI are orescrited when 
the entry is positively sfgned. Qnly CD ifs prescribed when the entry is 
unsigned. Only Cl fs prescribed when the entry is negatively signed. The CI entry 
on the first continuation card may be used when both CD and CI exist and area 
different. 


Cl references a CID entry on a coordinate system definition tnput card. when CI is 
entered as Q, interpolation is in basic rectangular coordinates. when Cl is not 
entered, interpolation is in the coordinate system that positions the fieid vertex 
and edge points. In this case, all vertex and edce points must be defined in the 
Same coordinate system (CP entries on GRID and/or cGRIOD cards). See Part [, 
Section 3.1 and Table 2. 


The 4 intersolation point components, in spherical coordinate systems, or 
the 9 interpolation point components, in cylindrical coordinate systens, may be 
offset by a +360° by appending a +S or a +€ to grid point entries. 


If an entry for GD is presert in a TRIA field, grid potnt 0D replaces grid point C 
in interpolating grid point positions from side CA. Grid point D is usually 
defined by an EGRID card. MM must be blank. See Part I, Section 3.5. 


[f both EL(XY) and E2(XY) are entered, the positions of generated grid points are 
determined by cubic fnterpolation from edge XY. [f E2(XY) is blank, positions are 
determined by quadratic interpolation. [f both EL(XY) and €2(X¥) are olank, 
positions are determined by linear interpolation. Either none or one edce point, 
E(XY) may be specified for each edge of a RIA field. &1(CD) positions 3 grid 
point on edge CA of IRIA fields. See Part [, Section 3.7. 


Intermediate continuation cards cannot be omitted even if all entries are dSlank. - 
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Input Data Card CGEN Generate Element Connection Bulk Oata Cards, Form 1 
Description: Generate element connection cards for element types J{sted under Remark 1 


Format and £xample: 


a 
eo fam foe [oe fat | | | be 
a ao 
ee a ae 









—_ Cr ===> 


Field Contents 

TYPE Any of the mnemonics listed under Remark 1 

FEIO The absolute vaiue of FEIO is the element identification number, EID, of first 
element in element set (Intecer # 0). See Remark 4. 

P{0 ID of an element property card or °GEN card (Integer > 0 or biank). Must be unique 
with respect to all other property card [0s and PGEN I0s. See Remark $. 

FIO ID of associated grid point field defined on a GRIOG or GRIDU card 

DIR One of the entries L, M, or N (Default entry is L). Direction of elements in the 


field. See Table 4 of Part I, Section 7.1. 


TH Material prooerty orientation specification (Real or Blank; or Integer > 0 and less 
than 1,000,000). If Real or 3lank, specifies the material property orientation 
angle in degrees. If integer, the orfentation of the material x-axis is along the 
projection onto the plane of the element of the x-axis of the co-ordinate system 
specified by the integer value. 


ETOL ID of element, below which no element connection cards are generated (Integer > 9 
or blank) 

EIDH ID of element, above which no element connection cards are generated (Integer > 0 
or blank) 


faene,'C,T0 Thickness of corner soints of QUAD and TRIA fields when TYPE is TRIA3, TRIA6, CUAD4S 
or QUAOS (Real or Blank). See Remark 6. - 


Remarks: 1. Element ‘types: R00, PLATEL, CFTUBE, TUBE, TRIA3, TRAPRG, TRIARG, TRIAG, TRIM6, 
TRIAX6, QUADS, SHEAR, QUADS, HEXAL, HEXAZ2, HEXA8, HEXA20 and HEX20. The mnemonics 
HEXAS and HEXA20 refer to HEXA elements with 8 or 9 through 20 nodes, respectively. 


(Continued) 


mh, 





ee 





CGEN (Cont. ) 


Elements are tabulated dy class in Table 3, of Part I, Section 7.0. Line class 
elements may be generated in all field types. Trfia3 and Tria6 class elements may 
be cenerated in TRIA and QUAD fields. Quad4 and Qiad8 class elements may be 
generated in TRIA, QUAD and HEX fields. Only QUAD4 and QUACD8 element types are 
valid in TRIA fields. TRIAI and TRIA6 elements are generated, in this case, along 
the field edge corresponding to the element direction. Hex8 and Hex20 class 
elements may only be generated in HEX fields. | 


MSGMESH automatically generates PENTA elements in place of HEXA and HEX20 elements, 
when adjacent element edges would have the same grid point identifiers. TRIA6 
elements replace QUADS elements and TRIA3 elements replace QUADS elements, when 
adjacent element vertices would have the same arid noint identifiers. Otherwise, 
MSGMESH suporesses the generation of an element which would reference identical 
grid point identifiers. Elements are generated only wnen all vertex points 
exist. Deleted midside nodes suppress the generation fo TRIM6 and TRIAX6 elements. 


EIDs are incremented when the FEID entry {is positive and decremented when the FEID 
entry {is negative. 


The PID entry may de blank for those elements whose element connection card permit 
a blank property iD. 


Element thicknesses may de generated for TRIA3, TRIA6, QUAD4 and CUADS8 elements 
located in TRIA and QUAD fields. Interpolation from the thicknesses at the field 
vertices is based on field topologies with equal spacings of thicknesses for both 
GRIOG and GRIDU fields. Unequal spacings of thicknesses are defined for GRIDG 
fields specified with unequal spacings of grid soints. 
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[Input Oata Card MATI Material Property Definition, Form 1 


Descriotion: Defines the material oroperties for linear, temperature-independent, isotropic 
materials 


Format and =xamole: 





ae aii = 





meld Contents 

MID Material identification number (Integer > 0) 

= Young's modulus (Real or blank) 

G Shear modulus (Real or blank) 

NU Poisson's ratio (-1.0 < Real < 0.5 or blank ) 

RHG Mass density (Real) 

A Thermal expansion coefficient (Real) 

TREF Thermal expansion reference temoerature (eal) 

GE Structural element damoing coefficient (Real ) 

Sac, 935 Stress limits for tension, compression, and shear (Real). (Used only to comoute 


margins of safety in certain elements; they have no effect on the computational 
Procedures. ) 


MCSID Material Coordinate System identification number (Integer > 0 or blank) 
Remarks: 1. The material identification number must be unique for all MAT1, MAT2, MAT3 and MATS 
cards. 
2. MATL materials may be made temperature dependent by use of the MATTI card. 


Ss The mass density, RHG, will be used to automatically comoute mass for all 
structural elements. ; 


4, Weight density may Se used in field 5 if the value 1/g fs entered on the PARAM card 
WIMASS, where g {ts the acceleration of gravity (see Section 3.1.5). 


5. MCSID must be nonzero if the CURVY module is used to calculate stresses or strains 
at grid points on plate and shell elements only. 


6. To obtain the damping coefficient, GE, multiply the critical damping ratio Clg, dy 
ZF 


(Continued } 
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10. 


MATL (Cont. ) 


Either — or G must be specified ({.e., nondlank). 

If any one of £, G, or NU is blank, ft will de computed to satisfy the identity 
FE = 2(1+NU)G; otherwise, values supplied by the user will be used. This 
calculation is only made for tnitial values of €, G, and NU. 

[f E and NU or G and NU are both blank, they will Doth be given the value 0.0. 


Implausible data on oe or more MATL cards will result in a warning messaqe. 
imolausible data is defined as any of €E ¢ 0.0, G ¢ 0.0, 0.5 < NU ¢ 0.0, or 


[l - 5 — c | > 0.01 (except for cases covered Sy Remark 9). 
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[npee Data Card MATSI Material Stress Dependence 


Descriotion: Specifies table references and material properties which are stress-denendent for 
use In material nonlinearity applications. This card will be activated {f MAT] with the same ‘ID 
1s being used in the material nonlinear solution sequence (66). 


Format and Examole: 


9 10 


vats ce ae = | 





MATS] 28 PEAS ile : lsagc 


7 


a ——— — | 


Field Contents 
MID Identification number of a MATL card (Integer > 0) 
Ri Identification number of a TABLES] card (Integer > 0) 
. VrYPt Type of material nonlinearity (8CD) NLELAST (Nonlinear elastic) or PLASTIC 


(Plasticity theory) 


| Work hardening slope (slope of stress vs. plastic strain) in units of stress 
(Real). For elastic-perfectly plastic cases, 8 2 0.0. For more than a single 
Slope in the plastic range, the stress-strain data must be supplied in a TABLES] 
card referenced in field 3. See Remarks 1 and 2. 


YF Yield functton (Integer); one of the following values: (Default = 1) 
l (von Mises) 
2 (Tresca) 
3 (Mohr-Coulomb ) 
4 (Drucker-Prager ) 
HR Hardening Rule (Integer); one of the following values: (Default = 1) 
l (Isotropic) 
2 (Kinematic) 
3 (Combined isotropic and kinematic hardening) 
LIMITI Parameter used in the yteld function specification (Real). See Remark 4. 
LIMIT2 Parameter used in the yield function specification for yield functions 3 and 4 


(Real). See Remark 4. 


8 3 








Remarks: Il. 





MATS1 (Cont.) 


If type = NLELAST, the stress-strain data given {n the TABLES] card will be use 
determine the stress for a given value of strain. The values 8, YF, HR, LIM 
and LIMIT2 will not be used in this case. ; 


[If type »* PLASTIC, either the table identification, R1, or the work harder 
slope, 8, may be gtven but not both. If the table [D0 is omitted, the 1 
hardening slope, 8, specified in field 3 is defined as 


; jee ere 
E E 


where £ is the elastic modulus and f’ {s the slope of the untaxial stress-str 
curve im the plastic region. 


If RL is given, TABLES! entries of stress-strain data must conform to the follow 
rules: 


as If type = NLELAST, the curve must pass through origin (x,20., y;70., for s 
value of 1). This option is valid only for line elements. 


b. If type = PLASTIC, only data for the first quadrant need be supplied. Th 
are two ootions available ta specify the data. In case 1 (defined as F@R!I 
tn the TABLES] card) the first point is at the origin, the second point is 
the yteld stress given in field 8 of the MATS1 card. The slope of the | 
jotning the origin to the yield stress must be equal to the value of € gt! 
on the MATL Bulk Data card. In case 2 (defined as FORM21 jin the TABLI 
card) the slopes of the curve in the plastic range are supplied. See ! 
TABLES] card description for additional details. 


LIMIT] and LIMIT2 are parameters used in the yield function definition as follows 


Yield Function LIMIT 


von Mises (1) Yield stress in tension, oa 
or Tresca (2) f 










LIMIT2 






Not used 








Mohr-Coulomd (3) 
or Orucker-Prager (4) 





Cohesion,c ({n stress units) Angle of internal 


friction » (in degree: 





See Section 2.14.1.2 {n the MSC/NASTRAN Application Manual for a detail 
description of yteld function and hardening rules. 
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MATSL (Cont.) 











NLELAST 
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[nout Oata Card TABLES) Matertal Property Table, Form | 


Description: Defines a tabular function for stress-dependent material properties such as the 
stress-strain curve and creep parameters 


Format and Examoie: 


] 7 3 
rs a SSSI 





Ee a OR a Ee 


Pata step spat t= — 









(etc. ) 
Field Contents 
[0 Table identification number (Integer > 0). 
[ Not used 
mrpety 4 Tabular entries (Real ) 


Remarks: 1. The xy must be in either ascending or descending order but not both. 
2. Jumps between two points (x, = Xi 1) are allowed, Dut not at the end joints. 
3. At least two entries must be present. 


4. Any x-y entry may be ignored dy placing the BCO string SKIP in aither of tne %wo 
fieids used for tnat entry. 


5: The end of the table is indicated by the existence of the 8CD string ENOT in aitner 
of the two fieids following the last entry. An error is detected if any 
continuation cards follow the card containing the end-of-table flag =NOT. 


6. TABLES] is used to input a curve in the form of 


Hs) Yt x) 


‘| 


aner2 X is input to the table and Y is returned. The table look-uo yr(x), x 2 X%, 
is performed using linear interpolation within the table and linear extrapolation 
Qutside the table using the last two end points at the appropriate table end. At 
jump points the average yr (x) is used. 


S6 





inout Data Card PL3ADZ Pressure Load on a ‘wo-Qimensional Structural Element 


Descrintion: Defines a uniform static pressure load applied to two-dimensional elements.  Inly 
, ontar, TRIA3, or TRIAX6 elements may have a pressure load applied to them via this card. 


Format and Example: 


1 2 3 4 5 5 7 g q 10 
PLOADZ SIo 9 Er9 pide) CIC EID FID | ce 
PLJAD2 75 aan a 4 | 16 2 | | 


Alternate Form: 


pLgan2 =| =~ SID > eto | “THRU” | loz | << | >< = 


Field Contents 

SO Load set identification number (Integer > 9) 

p Pressure value (Real) 

FID 

o1| Element identification number (Integer > 0; EIDL < £102) 


Remarks: 1. EID must Se 2 or blank for omitted entries. 


2. Load sets must de selected itn the Case Control Deck (LJAD=SID) to be used dy 
MSC/NASTRAN. 


Se. At least one positive £ID must be oresent an each PLUAD2 card. 


4, [f the alternate form is used, al] elements 2101 through E102 must de *two- 
dimensional. 


Ss. The direction of the pressure {s computed according to the riaht-hand rule using 
the grid point sequence specified on the element card. Refer to the PLZAD card. 


6. All elements referenced must exist. 


ee Continuation cards are not allowed. 
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V63 


[Input Data Card NLPARM Parameters for Nonlinear Analysis Control 
Description: Defines a set of parameters for nonlinear analysis iteration strategy 
Format and Example: 


l 2 3 4 5 6 7 8 9 10 
foram [8 [we | | erm] are [ware] ow | oon] 
ae ee 











Field Contents 

ID [dentification Number (Integer > 0) 

INC Number of fncrements (Integer > 0) 

OT Incremental time interval for creep analysis (Real > 0) (Default 2 9.; no creep) 
KMETH@O Method for controlling tangent stiffness updates (8CO 2 AUTO, ITER,SEMI ,LSON,AUTOON, 


or SEMION) (No default). See Remark 5. 


KSTEP Number of iterations before a matrix update. Used only for ITER and LSON 
methods. (Integer > 0 or blank) (Default = 1) 


MAXITER Limit on total iterations for each load increment (Integer * 9) (Oefault = 20) 

CONV Flags to select convergence tests (8CD = U, ?, W, or any combination of the 
letters) (Default 2 W). See Remarks 3 and 4. 

INTOUT Intermediate output flag (8CD = YES or NO) (Default 2 YES). See Remark §. 

EPSU,EPSP, Error tolerance limits for determining convergence {Real > 0.0) (Defaults 2 1.0&-2, 

ePSW 1.0&-3, 1.0£&-4, respectively) 

MAXDIV Limit on diverging iterations (Integer > 90) (Default 2 1) 

MAXQN Maximum number of Qudasi-Newton correction vectors to be saved on the data 5ase 


(Integer > 0) (Default 2 10) 


MAXLS Maximum number of linear search operations per vector iteration (Integer > 0) 
(Default 2 8) 


FSTRESS Fraction of effective stress (a) used ta determine the subincrement size in tne 
matertal routines (0. < Real < 1.0) (Default = 0.15). See Remark 9. 


(Continued) 
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Remarks: 


NLPARM (Cont. ) 


The NLPARM 3ulk Data card is requested by the Case Control card NLPARM = [D. Each 
solution subcase requires a loading condition and an NLPARM request. 


[In case of static analysis (DT = blank), INC is the number of equal subdivisions of 
the net load cnange defined for the SUBCASE. Applied toads, gravity loads, 
temperature sets, enforced displacements, etc., define the new loading condition. 
The differences from the previous case are divided by [NC to define the incremental! 
values. In case of creep analysis (DT > 0O.), INC is the number of time step 
increments. 


The test flags (U,P,N) and the error limits (EPSU, etc.) define the conversence 
criteria (U = Ofsolacement error test, 7? = Load equiltbrium error test, W = work 
error test). If the internally calculated error fractions are less than the 
requested limits, the iteration stoos, the results are srocessed, and the orogram 
Continues to the next load increment or sudcase. All requested error criteria must 
de satisfied for a “convergent solution”. 


[f convergence does not occur for a particular load increment, the number of 
iterations is limited to MAXITER. [If MAXITER is a positive number, the results 
will be processed and the program will continue with the next load increment. if 
MAXITER iS @ negative value, the solution will terminate if convergence is not 
acnieved. 


The basic nonlinear solution methods are the “QFGS" Quasi-Newton and -x70d1 fied 
“Newton-Raphson” processes wheredy the out-of-baiance nonlinear loads are neasured 
and recycled to solve for improved displacements. The convergence of this 
iteration process is generally improved by either performing “line searcnes” or dy 
frequent updates to the “tangent stiffness matrix" to account for current plastic 
and geometric effects. However, the matrix updates may de relatively costly and 
may be detrimental if the astimates are far from the correct answer. Thus, the 
user 1$ given the following KMETHODO options to provide control over a large range 
of problem types. 


a. [f the AUT2 option is selected, the program will automatically select tne most 
efficient strateqy Dased upon convergence rates. At each step the number of 
Iterations required to converge is estimated. [f this exceeds the MAXITER 
limit, or if the CPU time exceeds the time to perform a stiffness matrix 
update, the module will perform a stiffness update. If divergence occurs for 
MAXOLIY successive iterations, a stiffness update is also performed. 


Dis [If the SEMI option is selected, for each load increment the program will: (1) 
perform a single displacement iteration based upon tne new load; (2) odate 
tne stiffness matrix for the estimated geometry and material properties; and 
(3) resume the normal AUTD iteration method. 


Ga If the ITER option is selected, the stiffness matrix will be updated only 
after eacn KSTEP iteration. The count is reset to zero for each new matrix or 
load step. 


dis [f one of the AUTOQN, SEMIQN, or the LSON options jis selected tne “Quasi- 
Newton” and “Line Search” procedures are used. 3asically tne line searcn (LS) 
will scale the displacement increment vector to minimize the error at 2acn 
iteration. If the line search provides a substantial improvement, the Quasi- 
Newton (QN) operations use the local tangent of the error to correct the 
tangent matrix solution for subsequent iterations. Otherwise the operations 
behave much like the AUT@, SEMI, and ITER options. 


(Continued) 


89 








NLPARM (Cont.) 


Qutput requests for ELFQRCE and STRESS made in Case Control will be processed if 
(NTQUT = YES (the default value). [f (NTOUT = NQ, these output requests are 
Qrocessed for only the last load factor in a sudcase. Requests for OL[SP printout 
and structure plots are also processed during the nonlinear iteratton if PARAM, 
NLOISP, Ll is provided in the 8ulk Oata deck. 


[f a diverging solution {s encountered for MAXOIVY successful {terations, the AUTG, 
AUTGQN, SEMI, and SEMIQN options will perform a tangent stiffness update. La 
after the stiffness update, the solution is stil! diverging, the program 
terminates. For {TER and LSQN options, diverging cases cause immediate 
termination. 


The unit of the incremental time interval under OT field must de consistent with 
the unit used on the CREEP card to define the creeo characteristics. Total creep 
time for the sudDcase may be found from OT multiplied by the value in the [NC. 


The number of subincrements in the material routines (elastonplastic and creep) are 
determined such that the subincrement size is approximately FSTRESS * 3. 


FSTRESS is also used to establish a tolerance for elastoplastic material in the 
matertal routine, {i.e., 


6 in yteld function ¢ FSTRESS * yteld stress : 


[f the limit is exceeded at the converging state, the program wil! exit with a 
fatal error message. Qtherwise the stress state is adjusted to the current yield 
surface, resulting in 5 = Q. 
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Input Data Card SPCG Constraint Set Generator, Standard Form 


Description: Generates SPC1 cards that constrain a subset of grid points in a grid point field 


Format and Examole: 





SID ID number of single-point constraint set (Integer > 0) 

FID ID number of a GRIDG or GRIDU field 

c Component numbers (any combination of the digits 1-6 with no imbedded blanks) in 
the displacement coordinate system (CD entry on GRIDG card) 

01,02 Location number pair of order-independent pivot points that define a subset of grid 
points in FID. See Part I, Section 2.5. 


* 


Remarks: l. See Part I, Section 5.0 for a discussi{on on the SPCG input card. 


2. SPC entries on generated SPCl cards are made only for active grid points. See 
Part I, Section 10.1. 


elee 





[Input Data Card ODLGAD Dynamic Load Combination (Superposition) 


Nescription: Defines a dynamic loading condition for frequency response or transient response 
problems as a linear combination of load sets defined via RLJADL or RLGAD2 cards (for frequency 
response) or TLJADL or TLYAD2 cards (for transient response) 


Format and Examoles: 


1 2 3 4 ) 5 7 i 9 10 
SE RS aa a A a 
ng | oi | iio | zof 6 | -2of 7 | eof is [a | 
ee ae ee a A a A 
LED ees ge eee a ee a 












-etc.- 
Field Contents 
S10 Load set identification number (Integer > D) 
5 Scale factor (Real) 
Bs Scale Factors (Real) 
= Load set identfication numbers defined via card types enumerated above 


(Integer > D) 


Remarks: 1. The load vector being defined by this card is given by 
iPim= Ss) S{P.} 
i 


cs The Li must be unique. 
ie SID must be unique from all Li. 


4, Nonlinear transient loads may not be included; they are selected separately in the 
Case Control! Deck. 


5. Linear load sets must be selected in the Case Control Deck (DLIJAD=SID) to be used 
by MSC/NASTRAN. 


6. A DLJAD card may not reference a set identification number defined by another OLGAD 
card. ’ 


7.  TLGADL and TLGAD2 loads may be combined only through the use of the OLJAD card. 
8. RLYADL and RLZAD2 loads may be combined only through the use of the OLGAD card. 


9. SID must be unique for all TLGAD1L, TLGADZ, RLGADL and RLGADZ cards. 


BZ 








Input Data Card TLGAD1 Transient Response Dynamic Load, Form 1 


Descriotion: Oefines a time dependent dynamic load or enforced motion of the form 


(OG) ea Ame (toa |} 


for use in a transfent response oroblem 


Format and Examole: 


7 2 3 4 5 6 
for [oe |. [| [| «Pee 
roo [s [7 |] =] sl 1 






\ 


| 
X 
I 








Field Contents 

SID Set identification number ([nteger > 0) 

L Identification number of DAREA card set or a static load card set or a thermal load 
set (in heat transfer analysis) which defines A (Inteser > 0) 

M Identification number of DELAY card set which defines <= (Integer > 0) 

TYPE Defines the nature of the dynamic excitation (Integer 0, 1, 2, 32 or dBlank). (See 
Remark 2.) 

TIO Identification number of TABLED! card whicn gives Flt - t}) (Integer > 9) 


Remarks: 1. If M fs zero, twill be zero. 


2. The nature of the dynamic excitation {s defined in field S in accordance with the 
following tabdle. 


0 or blank 








Force or Moment 





Enforced Displacement 
Enforced Velocity 


Enforced Acceleration 


see Section 2.7 of the MSC/NASTRAN Application Manual regarding tne use of the 
enforced motion options. Note that “large masses” must be used for enforced 
motion. For heat transfer problems, Field 5 must be Ddlank. 


(Continued) 
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TLOADL (Conc. ) 


Dynamic load sets must be selected {n the Case Control Deck (OLJAD=SID) to be used 
by MSC/NASTRAN. 


TLYAD1 loads may de combined witn TLIADZ loads only dy specification on a DLGAD 
card. That is, the SID on a TL@ADL card may not be the same as tnat on a TLIADZ 
card. 


S(D must be unique for all TLJAOL, TL@AD2, RLGJAOL and ALJAD? cards. 


Field 3 may reference sets containing QHBDY, QBDY1, QBDY2, QVECT, and QV@L cards 
wnen using the heat transfer option. 


[f the heat transfer option is used, the referenced QVECT data card may also 
contain references to functions of time, and therefore A may de a function of time. 


Fourter analysis will be used if this {s selected in an aeroelastic response 
9robdlem. 
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Inout Oata Card LSEQ Static Load Set Definition 


Descriotion: Oefines a seouence of static load sets to be applied to the structur3l model]. The 
Toad sets may be referenced by dynamic load cards. 


Format and Examole: 


1 2 3 4 5 6 7 8 9 10 
fe [| pea] om] m) | | | | 
S| to | mo | wo] o |. ||. |. | 





Field Contents 

$10 Set identification of the set of LSEQ cards (Integer > 0} 

DAREA The DAREA set Identification assigned to this static load vector (Integer > 0) 

LID Load set identification number of a set of static load cards (Any card that may be 


referenced dy the LZAD Case Control card) (Integer > 0 or blank) 

p80) Set identification of a thermal load set (Any card that may de referenced dy the 
TEMP(LBA0) Case Control card) (Integer > 0 or blank) 

Remarks: 1. The above cards will not be used unless selected in the Case Control Deck with a 


L@AOSET card. 


2. This card is available only in superelements and in Solution Sequences 26, 27, 30 
andes . 


ar The number of static load vectors created for each suoerelement is the number of 
unique DAREA [Ds on all LSEQ cards in the bulk data. 


4, The DAREA identification assigned to the static load vectors may be referenced dy 
RLZAODL, RLZADZ, TL3AD1 and TLOAD2 cards. 


5S. Element data recovery fer thermal loads is not currently imolemented in dynamics. 


6. DAREA set identification numbers should be unique with respec 72 311] static load 
set identification numbers. 
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[noput Data Card 


Description: 
dependent dynamic loads 


TABLEOL Oynamic Load Tabular Function, Form 1 


Defines a tabular function for use in generating frequency-dependent and time- 


Format and Examole: 


4 5 6 7 8 9 19 
<P te i A ae 
TABLED L | | “ABC 
— a Te | 


ape eee 





X4 oF 4 


Remarks: 


Soo etc.- 


Contents 
Table identification number (Integer > 9) 


Tabular entries (Real ) 


The Xy must de tn etther ascending or descending order but not Sotn. 
Jumps between two points (x, = X41) are allowed, but not at the end points. 
At least two entries must Se present. 


Any x-y entry may be ignored dy alacing tne 3CD string “SKIP” in aitner of “Ne two 
fieids used for that entry. 


The end of the table {s itndicated by the existence of tne 3C0 string “SNOT" in 
either of the two fields following the last entry. Am error is detected if any 
continuation cards follow the card containing the end-of-table flag “ENOT”. 

Eacn TABLED mnemonic infers the use of a specific aigoritnm. For TABLED] type 
tables, this algorithm is 


Y = y,(X) 


where X is input to the table and Y is returned. The table look-up yy. x}, x = X, 
is performed using linear interpolation within tne table and linear extrapolation 
outside the table using the last two end points at the approoriate table end. At 
jump points the average yr(x ) is used. There are no error returns from this taole 
look-up procedure. = 


Linear extrapolation {s not used for Fourier Transform methods. The function is 
zero outside the range. 
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Inout Data Card 


Descriotion: 


= 


TSTEPNt Timm Step and Nonlinear Solusion Tontzry 


Provices Parametric controls and cata for Nonlinear Transient Anelysis 


rPormat and exzamole: 





] 
Ser 


om | 
= 
a 


TSTEPNL 





MAXITEP 


CONG 


2PSU cP SP EPS 


MAXD (4 


MAXON 


NOT ae W | wmeTagD | «sTep | wexttes | 
250 an =| ni | > | wie | | ay 
PSU EPSP epsi | masory | maxon | manus | #s72ss | 
01 D se Rood | ane Lenco | 

Cont ants 





[dentiftcarion 


Number of *imm steos 


Time increment (eal > 9.9} 

Skio factor for sstout (every MEN Step 011] Ge caved for osstout} (1 
mathed for controlling iterations and tanaent matris uscates (30) = 
(No default) 

Numer of steps before a matrix update (integer > 0) (Befawiz * 2) 
Limit om coral iterations for eacn tise increment (letecer > 9) (De 


to selecs convernence 


S 
ers or olan) 


AW Error tolerance iimits fear detecmioing <enversence [Peal » 2.93 
~ - # ~ 4 \ 
(Defaults = 1.06-2, 1.0€-2, 1.0€-4, cmspectiveiy 
Limit on the mumper of diwerzing iteragtoms 2f each time ang 2 
stees {Integer > 0) (Oafauit = 2) 
Maximum numcer cf Quesi-Newton correction vectors 19 we sawed gn 


({aweger > 9) (Befaylzt = 29) 


t ~~ ° 
Maximum nummer of ‘Line Searcn iterations for ‘LS0M aptior 
(Defayiz = 2) 
-F2 < ¢¢ oer ts > - ad . aa a 
Fpaction «f offecztiwe stress 7) £522 (3 AL ermine Ie wl) ads seer 
° Ba « fr 9 fe 4 = 
material routines (0. < Amal < 2.0) {Defawig, = 9.15). See 2emere 20 
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Remarks: 


Ls 


NASTRAN DATA OECK 


TSTEPNL (Cont. ) 


The TSTEPNL Bulk Oata card {ts selected by the Case Control card TSTEPNL = [0. Each 
sudcase (residual superelement solutions only) requires a TSTEPNL card and either 
applied loads vita TLQADi data or initial values from a previous subcase. Multiple 
subcases are assumed to occur sequentially in time with the initial values of time 
and displacement condittons of eacn subcase defined by the end conditions of tne 
Drevious subcase. 


OT, NOT, N@ are the time increments, 4t, the total numder of steps \,, and the out- 
gut interval, Nos respectively. The value of time, T,, at each increment n will be 


t, « Cy + nat (ae oa eorre Ns 


where ty is the initial value of time which corresponds to the last step from tne 
previous subdcase. 


For printing and plotting the solution, data will be outout at steps n = Q, Ny, 
2gy--eNe The Case Control card JTIME may also de used to Control the outsut 
points. 


The basic solution method {s the “Newmark 8eta” transient numerical method along 
with tne "Newton-Raphson" metnod for controlling nonlinear effects. Equilidrium 
solutions with respect to inertial, damping, elastic, and nonlinear loads are 
obtained at each time step much the same as in the statics case. 


The KMETHQO options provide the user with tnree optional sxonlinear iteration 
metnods. These are: 


AUT9 provides for automatic stiffness matrix updates to improve poor 
convergence. The KSTEP value is ignored. 


TSTEP provides for a matrix update at every xsTep®” increment of time 
(typically a large number). 


The stiffness matrix is always updated for a new subcase or restart. OInly one 
update {s allcwed for eacn time steo and the iteration process returns to the 
Qrevious time step after an update. 


MAXITER defines the number of static nonlinear load iterations (corrections) to de 
performed at each time step. At least two iterations are necessary *o cetermine 
the conversence of the displacements (see below). 


The test flags (U,P.Ww) and tne error limits (EPSU, etc.) define the canvergencea 
criteria (U 2 Displacement error test, P = Load equilibrium error test, W 2 Work 
error test). [f the internally calculated error fractions are less than the 
requested limits, the fteration stoos, the results are processed, and the program 
continues to the next time steo or subdcase. 


Se: 





NOfe 


TSTEPNL (Cont. } 


If no tests are requested or mo convergence occurs, the number of iterations is 
limited to MAXITER. The results will be processed and the program will continue 
with the next time step. 


MAXQN provides tne control over diverging solutions. At each time steo, the 
solution is allowed to diverge until, at the MAXDIVY iteration, a matrix update is 


requested. If the average convergence parameter i, indicates divergence for ‘“AXDLY 


t 
successive time steps, the transient solution stops and exits for data recovery. 


The number of subincrements in tne material routines (elastoplastic and creep) are 
determined such that the subincrement size is approximately FSTRESS * a. 


FSTRESS {1s also used to establish a tolerance for elastoplastic material in the 
material routine, ji.e., 


6 in yteld function < FSTRESS * yield stress : 


If the limit {is exceeded at the converging state, the program will EXIT with a 
fatal error message. Otherwise the stress state is adjusted to the current yield 
surface, resulting tn 6 = Q. 
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